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I .  INTRODUCTION 

The  objective  of  the  regional  discrimination  research 
project  is  to  systematically  assess  the  potential  of  regional 
phases  for  discriminating  earthquakes  and  nuclear  explosions 
with  particular  emphasis  on  defining  the  applicability  of  such 
discriminants  to  events  in  various  regions  of  the  U.S.S.R. 

The  research  program  includes  two  distinct  approaches.  The 
first  employs  more  traditional,  time  domain  amplitude  and 
period  measurements  to  compare  the  relative  excitation  of  dif¬ 
ferent  regional  phases  from  earthquakes  and  explosions.  The 
second  approach  focuses  on  the  evaluation  of  the  capabilities 
of  other  measurement  or  processing  techniques,  such  as  spectral 
ratios  or  narrow-band  filtering,  to  identify  seismic  sources 
using  regional  phases. 

The  principal  technical  phases  of  the  regional  discrimina¬ 
tion  research  effort  are  as  follows: 

•  Review  and  assess  the  time  and  frequency 
domain  characteristics  of  regional  seismic 
phases  and  evaluate  their  dependency  on 
source  type,  propagation  path  and  station 
environment . 

•  Develop  a  tentative  set  of  source  identi¬ 
fication  criteria  based  on  information 
from  regional  phases. 


i 

i 


» 


e  Test  proposed  regional  discriminants  on 
observed  data  to  determine  their  capability 
to  identify  earthquakes  and  explosions  in 
various  tectonic  and  geologic  settings. 

In  Section  II,  the  Technical  Discussion,  we  describe 
our  research  accomplishments  to  date  on  the  current  project. 

The  impact  of  our  current  results  on  research  efforts  antici¬ 
pated  under  the  continuation  contract  is  outlined  in  Section 
III.  A  description  of  the  results  of  our  analysis  of  seismo¬ 
grams  for  Eurasian  events  recorded  at  AEDS  stations  is  presented 
in  a  separate  appendix  to  this  report. 
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II.  TECHNICAL  DISCUSSION 


2.1  REVIEW  OF  LITERATURE  ON  REGIONAL  PHASE  PROPAGATION 

FROM  RUSSIAN  EVENTS 

2.1.1  Introduction 

Much  of  the  interest  in  the  possible  use  of  regional 
phases  in  discriminating  between  underground  nuclear  explo¬ 
sions  and  earthquakes  stems  from  observational  data  which 
suggests  that  in  certain  tectonic  environments  such  phases 
produce  the  largest  amplitude  out  to  between  20°  and  30° 
from  the  source.  For  smaller  events,  regional  phases  in 
this  distance  range  may  produce  the  only  detectable  signal. 

In  the  following,  we  summarize  the  state  of  our  understanding 
of  regional  phase  characteristics  in  the  U.S.S.R.  as  derived 
from  observational  data.  The  characteristics  include  propaga¬ 
tion  and  attenuation  of  regional  phases  as  well  as  studies 
involving  their  application  to  discrimination. 


2.1.2  Propagation  of  Regional  Phases  Crossing  the  U.S.S.R. 

Our  understanding  of  the  propagation  of  regional  phases 
in  the  interior  of  the  U.S.S.R.  has  been  hampered  by  the  lack 
of  in-country  seismic  data.  Although  some  information  on  re¬ 
gional  phases  from  earthquakes  in  Russia  is  now  available,  no 
similar  information  exists  for  underground  nuclear  explosions 
except  that  recorded  at  stations  outside  the  U.S.S.R.  How¬ 
ever,  the  usefulness  of  the  latter  data  is  itself  limited 
because  of  the  lack  of  correlative  earthquake  data  with  sources 
near  the  explosion  sites  and  our  limited  understanding  of  how 
regional  phase  propagation  is  affected  by  geologic  structure 
along  the  path,  particularly  for  paths  crossing  the  tectoni¬ 
cally  complex  regions  of  Russia's  southern  and  eastern  borders. 

In  this  description  we  will  be  concerned  with  the  four 
regional  phases  which  have  received  the  most  attention  in  the 
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literature:  P  ,  P  ,  L  and  to  a  lesser  extent  S„ .  Probably 
n  g  g  n 

the  most  abundant  type  of  information  available  on  these  phases 

is  velocity  of  propagation.  This  type  of  information  comes 

from  a  variety  of  sources  including  single  and  multiple  station 

observations  of  earthquakes  at  regional  distances  (cf.  Bath, 

1954,  1957,  1959;  Pec,  1961;  Nersesov  and  Rautian,  1964;  Gupta, 

et  al. ,  1980)  as  well  as  refraction  surveys  (.for  p^  and  P  ) 

with  large  spread  lengths  such  as  those  encountered  in  Deep 

Seismic  Sounding  studies  which  have  been  conducted  in  many 

areas  of  the  U.S.S.R.  (cf.  Volvovsky,  1973;  Piwinskii,  1979). 

Figure  1  shows  some  typical  regional  phase  travel  time  curves 

and  velocities  for  southern  Russia  (.Pamir-Lena  River)  developed 

by  Nersesov  and  Rautian  (1964)  from  earthquake  data.  Gupta, 

et  al.  (1980)  show  similar  regional  phase  travel  time  curves 

for  earthquakes  and  explosions  in  western  Russia  as  observed 

at  WWSSN  stations  surrounding  the  Soviet  Union  at  distances 

from  5°  ^o  30° .  Their  plot  for  shallow  earthquakes  is  shown  in 

Figure  2.  For  all  the  events  studied  the  data  indicate  Pn  or  P 

velocities  varying  with  distance  from  8.26  to  8.49  km/sec  for 

the  distance  range  from  5°  to  16°,  from  9.8  to  10.38  km/sec  for 

the  distance  range  from  16°  to  25°  and  from  12.19  to  12.96 

km/sec  for  the  distance  range  beyond  25°.  Velocities  for  Pg 

were  found  to  be  between  5.52  and  5.66  km/sec,  for  S  between 

n 

4.67  and  4.71  kra/sec  and  for  Lg  from  3.49  to  3.57  km/sec. 

These  velocities  are  in  good  agreement  with  those  shown  in 
Figure  1  for  southern  Russia  and  appear  to  be  representative 
of  regional  phase  velocities  found  throughout  most  of  the 
Soviet  Union. 

More  important  to  the  discrimination  problem  are  the 

amplitudes  and  attenuation  of  the  regional  phases.  The  most 

abundant  published  data  on  regional  phase  amplitudes  in  the 

U.S.S.R.  appears  to  be  that  regarding  L^;  some  information  is 

available  on  P  and  P„  but  very  little  dealing  with  S  . 

n  g  n 

Piwinskii  and  Springer  (1978)  summarized  the  observations  of 
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several  authors  (viz.  Bath,  1954,  1957,  1959;  Payo,  1960;  Utsu, 
1958,  1960;  Pec,  1961,  1962;  Savarensky  and  Valdner,  1960; 
Sikharulidze,  1963;  Ruzaikin,  et  al.,  1977;  Pshenikov,  1961; 
Valdner  and  Savarensky,  1961;  and  Rijikova,  1966)  regarding 
the  propagation  of  phases  across  Europe  and  Asia.  Although 
the  data  reported  are  for  the  most  part  qualitative,  they  do 
provide  a  general  framework  for  understanding  the  efficiency 
of  Lg  propagation  in  these  areas.  Taken  as  a  whole,  the  data 
suggest  that  propagates  efficiently  throughout  the  platform 
areas  and  older  fold  systems  of  Asia  and  Europe  (cf.  Figure  3). 
The  younger  fold  belts  of  the  Alpine-Himalayan  system  and  re¬ 
gions  of  oceanic  crust  produce  erratic  propagation  or  act  to 
completely  block  transmission. 

The  study  by  Ruzaikin,  et  al^.  (1977)  appears  to  have 
particular  significance  to  assessing  the  potential  value  of  Lg 
from  events  in  south-central  Russia  recorded  at  regional  dis¬ 
tances  to  the  south.  Figure  4  shows  their  qualitative  assess¬ 
ment  of  Lg  from  events  throughout  southern  Asia  recorded  at 
two  stations  in  southern  Russia.  Darker  circles  indicate  that 
the  earthquakes  at  those  epicenters  produced  weak  or  no  Lg  at 
the  recording  stations  shown  as  triangles.  The  poor  Lg  propa¬ 
gation  from  events  to  the  south  is  attributed  to  strong  attenua¬ 
tion  or  interruption  of  the  wave  guide  in  the  Himalayan  fold 
system  or  Tibetan  Plateau  (cf.  Figure  3).  Using  reciprocity 
it  may  be  possible  to  infer  that  events  from  the  Soviet  nuclear 
testing  area  in  east  Kazakh  would  be  expected  to  produce,  at 
best,  weak  Lg  phases  at  recording  stations  near  the  darker 
circles . 

Studies  by  Gupta,  et  al.  (1980)  were  also  in  agreement 
with  these  findings.  They  reported  on  the  efficiency  of  Lg 
propagation  in  southern  Asia  to  stations  at  Kabul,  Afghanistan 
and  Meshed,  Iran.  Their  findings  are  summarized  in  Figure  5. 

Clear  L  phases  were  found  for  nearly  all  paths  from  events 
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north  and  northeast  of  the  stations  to  ranges  of  30°.  However, 
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Figure  3.  Tectonic  subdivision  of  Asia. 


Figure  4.  Efficiency  of  Lg  propagation  in  the  vicinity  of 
Russia's  southern  border  region.  Earthquake 
epicenters  indicated  by  circles;  stations  by 
triangles  (after  Ruzaikin,  et  al. ,  1977). 

8 


SYSTEMS.  SCIENCE  AND  SOETWAHE 


.*  0  - 


for  paths  crossing  the  southern  Caspian  Sea  and  along  the 
Himalayan-Alpine  fold  belt,  coincident  with  Russia's  southern 
border,  behavior  of  the  Lg  phase  was  more  erratic  with  fre¬ 
quent  cases  of  weak  or  no  phase  recorded. 

Similar  qualitative  information  on  the  propagation 
efficiency  of  other  regional  phases  is  not  generally  available. 
However,  considering  that  Lg  and  Pg  phases  are  commonly  asso¬ 
ciated  with  the  presence  of  a  granitic  layer  in  continental 
type  crust,  the  factors  affecting  Lg  attenuation  or  interrup¬ 
tion  of  the  wave  guide  might  be  expected  to  similarly  affect 
Pg.  Although  thickening  of  the  crust  under  the  Himalayas 

could  also  affect  P  and  S  phases,  the  attenuation  would  not 

n  n 

necessarily  be  expected  to  be  as  severe  as  that  for  Lg  and  Pg. 

More  specific  information  on  regional  phase  amplitudes 

and  attenuation  from  events  recorded  at  seismic  stations  within- 

the  U.S.S.R.  has  been  published  for  some  areas.  Figure  6  shows 

amplitude  attenuation  curves  for  L  ,  P  ,  P  and  P  phases  ob- 

g  g  n 

served  along  a  line  of  stations  in  southern  Russia  from  earth¬ 
quakes  in  the  Pamir  Mountains  (top)  and  Baikal  area  (bottom) 

(cf .  Nersesov  and  Rautian,  1964) .  In  general,  P_  and  L  are 

g  g 

seen  to  exhibit  similar  attenuation  behavior  with  distance. 

The  Lg  amplitudes  larger  than  those  of  any  of  the  P  phases 
out  to  a  distance  range  of  2500  km  where  they  are  more  severely 
attenuated.  It  is  also  noteworthy  that  Pg  data  are  not  re¬ 
ported  beyond  1200  km;  at  that  distance  their  amplitudes  become 
about  the  same  level  as  the  other  P  phases.  The  attenuation  of 
Pn  phases  particularly  at  shorter  distances  appear  more  erratic. 

Nuttli  (1980,  1981)  has  recently  published  more  detailed 

data  on  L  attenuation  from  Russian  events  recorded  at  seismic 
9 

stations  outside  the  U.S.S.R.  In  his  analysis,  Nuttli  has  cal¬ 
culated  an  effective  attenuation  rate  expressed  in  terms  of  a 
specific  quality  factor  QQ  for  several  nuclear  explosions  and 
some  earthquakes  in  Russia  observed  along  paths  to  the  northeast 
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and  south.  Although  some  questions  still  exist  regarding  the 
reliability  of  the  single  station  techniques  used  by  Nuttli  to 
estimate  attenuation  in  this  study,  the  results  are  in  general 
agreement  with  previously  cited  findings.  These  results  are 
indicated  in  Figures  7  through  10.  Referring  again  to  Figure 
3,  the  data  indicate  high  Q  (low  attenuation)  across  the  plat¬ 
form  regions  but  low  Q  (high  attenuation)  over  most  paths  to 
the  south  into  Iran  and  Afghanistan.  However,  some  of  the 
paths  to  the  south  crossing  the  Himalayan-Alpine  fold  belt 
appear  to  have  higher  than  expected  Q  estimates.  This  is 
particularly  notable  for  the  east  Kazakh  explosions  shown  in 
Figure  10.  Here  the  Qq  estimates  suggest  attenuation  levels 
intermediate  between  those  of  the  eastern  and  western  United 
States.  For  purposes  of  comparison,  the  QQ  for  one  hertz  Lg 
waves  is  estimated  as  1200  for  the  central  United  States, 

900  for  the  eastern  United  States  and  200  for  the  western 
United  States.  The  findings  of  Figure  10  do  not  seem  to  be 
in  agreement  with  the  results  of  Ruzaikin,  et  al.  (1977) 
shown  above  in  Figure  4  which  suggested  nearly  complete 
blockage  of  along  some  of  these  paths. 


Discrimination  Using  Regional  Phases 


With  regard  to  their  use  in  discrimination,  much  of 

the  research  on  regional  phases  has  focused  on  measurements 

of  the  relative  excitation  of  the  P  and  L  phases.  The 

g  g 

possibility  of  differences  in  excitation  level  for  earthquakes 
and  explosions  has  evolved  primarily  from  observations  in  the 
eastern  United  States.  Blandford  and  Hartenberger  (1978) 
suggested  that  the  ratio  of  the  maximum  amplitude  before  Sn 
to  the  maximum  amplitude  after  Sn,  as  measured  on  short  period 
vertical  instruments,  was  a  good  discriminant  between  earth¬ 
quakes  and  explosions.  Subsequent  studies  (e.g.,  Pomeroy, 
1980?  Gupta,  et  al. ,  1980;  Gupta  and  Burnetti,  1981)  have 
attempted  to  apply  this  technique  to  Russian  events  recorded 
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Figure  7 


Attenuation  of  L„  phases  in  terms  of  Qg  for 
paths  in  western  and  southern  Asia  from 
earthquakes  in  the  eastern  Caucasus.  Event 
epicenters  indicated  by  (Q) ;  station  loca¬ 
tions  by  small  circles. 
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8.  Attenuation  of  L„  phases  in  terms  of  Qg  for 
paths  in  central’and  southern  Asia  from  events 
in  west  Kazakh.  Location  of  explosions  indi¬ 
cated  by  (X) ;  earthquake  locations  by  (Q) ; 
station  locations  by  small  circles. 
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Figure  9.  Attenuation  of  Lg  phases  in  terms  of  Q0  for 

paths  in  Central  Asia  and  Northern  Europe  from 
explosions  at  Novaya  Zemlya.  Event  epicenters 
indicated  by  (x) ;  station  locations  by  small 
circles  (after  Nuttli,  1980) . 
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Figure  10.  Attenuation  of  Lg  phases  in  terms  of  Qq  for 
paths  in  Central  and  Southern  Asia  from  ex¬ 
plosions  in  east  Kazakh.  Event  epicenters 
indicated  by  (x) ;  station  locations  by  small 
circles  (after  Nuttli,  1980) . 
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outside  the  country.  In  these  cases,  the  maximum  vertical 
amplitude  after  Sn  has  generally  been  interpreted  to  corre¬ 
spond  to  Lg  while  that  prior  to  Sn  is  usually  attributed  to 
Pg  as  suggested  by  Figure  6  above.  Pomeroy  (1980)  shows 
peak  Lg  versus  maximum  P  wave  amplitudes  for  earthquakes  and 
explosions  in  central  and  western  Russia  (cf.  Figure  11). 

In  the  central  U.S.S.R.  the  Lg  to  P  wave  amplitude  ratios 
are  generally  less  than  one  for  explosions  and  greater  than 
one  for  earthquakes.  Although  there  is  considerable  scatter 
in  the  explosion  data,  the  results  appear  to  offer  some  hope 
for  use  of  such  measurements  in  discrimination.  The  results 
for  the  western  U.S.S.R.  are  not  so  promising.  Both  the 
earthquake  and  explosion  data  show  a  high  degree  of  scatter 
and  there  does  not  appear  to  be  much  consistency  with  regard 
to  Lg  to  P  amplitude  ratios  being  less  than  or  greater  than 
one.  Pomeroy  has  not  identified  propagation  paths  for  in¬ 
dividual  data  points  so  it  is  unclear  to  what  degree  attenua¬ 
tion  differences  might  be  contributing  to  the  observed  scatter. 
Gupta  and  Burnetti  (1981)  compared  the  maximum  amplitude  prior 
to  Sn  (maximum  P)  and  maximum  amplitude  after  Sn  (generally 

attributed  to  L  )  and  also  the  maximum  amplitude  in  the  P 
g  n 

window  to  the  maximum  amplitude  in  the  Lg  window  for  events 
in  the  southern  U.S.S.R.  recorded  at  Kabul  and  Meshed  (cf. 
Figure  12) .  Although  they  found  some  separation  between  the 
mean  trends  of  the  data  for  the  earthquakes  and  explosions, 
they  show  considerable  scatter  for  individual  data  points. 

It  seems  clear  from  the  available  peak  amplitude  data,  such 
as  that  described  here,  that  some  method  for  reducing  the 
scatter,  possibly  by  corrections  for  attenuation  or  source 
effects,  is  needed  if  such  measurements  are  to  evolve  into 
a  reliable  discrimination  tool. 

Pomeroy  (1979,  1980) ,  based  on  consideration  of  the 
excitation  characteristics  of  higher  mode  surface  waves, 
hypothesized  that  the  peak  of  the  Lg  motion  could  be  expected 
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to  have  a  lower  group  velocity  for  shallow  events  than  for 
deep  events.  Thus,  assuming  that  explosions  are  generally 
somewhat  shallower  than  earthquakes,  the  group  velocity 
corresponding  to  the  peak  amplitude  in  explosions  would 
be  expected  to  be  lower  than  that  for  earthquakes  from  the 
same  area.  Pomeroy  has  tested  this  hypothesis  using  two 
approaches:  one  based  on  a  single  measurement  of  the  veloc¬ 
ity  corresponding  to  the  peak  motion  and  the  other  based 
on  a  computation  of  energy  in  a  high  and  low  velocity  window 
of  the  Lg  motion.  Pomeroy's  results  for  the  eastern  United 
States  appear  to  agree  with  the  proposed  hypothesis;  however, 
the  application  to  the  U.S.S.R.  is  less  convincing.  Figure 
13  shows  the  group  velocities  determined  from  peak  motions 
for  events  in  the  U.S.S.R.  The  relatively  small  earthquake 
sample  is  seen  to  have  group  velocity  values  both  above  and 
below  the  bulk  of  the  explosion  data.  Additional  refinements 
using  data  with  common  source  regions  and  transmission  paths 
show  somewhat  more  optimistic  results;  but  more  complete  test¬ 
ing  using  a  larger  earthquake  sample  is  needed  to  assess  this 
technique. 

2.1.4  Summary 

Several  observations  can  be  made  regarding  regional 
phases  from  Russian  events: 

•  Propagation  velocities  associated  with  regional 
phases  in  the  Soviet  Union  are  about  7.9  to  8.3 
km/sec  for  Pn  (somewhat  higher  at  larger  dis¬ 
tances)  ,  5.6  to  6.0  km/sec  for  Pg,  about  4.7 
km/sec  for  Sn  and  about  3.5  km/sec  for  Lg. 

•  Regional  phases,  and  in  particular  Lg,  appear 
to  propagate  efficiently  along  most  paths 
within  the  Soviet  Union  -  i.e.,  platform  areas 
and  older  fold  systems. 


20 


SK5ffWS,  SCtCNCC  ANO  SOrTWAHE 


•  • 


•  • 


% 


06 

W 

(fl 

D 
O 
(N 


•  ■ 


O 

m 


(3SS/UIJ()  AXIOOTJA  dnOHD  (Z)  1 


21 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


* 


DISTANCE  (deg) 

Figure  13.  Group  velocities  corresponding  to  peak  Lq  amplitudes  for 
explosions  (circles)  and  earthquakes  (squares) . 


•  The  regional  phases  Lg,  and  probably  P~,  do 
not  propagate  efficiently  across  the  rlgion 
to  the  south  of  Russia  —  i.e.,  the  younger 
fold  system  of  the  Alpine-Himalayan  orogeny. 

•  Data  from  at  least  some  areas  of  the  U.S.S.R. 
indicate  that  L„  produces  signal  amplitudes 
from  earthquakes  larger  than  any  P  phases  out 
to  distances  of  2500  km  or  so.  P_  phases  from 
the  same  area  are  not  reported  beyond  about 
1200  km  where  their  amplitude  level  is  equiva¬ 
lent  to  that  of  other  P  phases. 

•  The  effective  attenuation  for  one  hert2  Lg 
waves  in  the  platform  and  older  fold  regions 
of  central  Russia  is  low,  approximately 
equivalent  to  that  in  the  eastern  United 
States.  Attenuation  in  the  southern  border 
region  is  high,  more  like  values  measured 

in  the  western  United  States. 

•  Discriminants  based  on  measures  of  peak 
amplitudes  of  Lg,  Pg  and  Pn  show  some 
evidence  of  separation  between  the  means 
of  earthquakes  and  explosion  sets,  but 
the  scatter  in  the  data  does  not  currently 
permit  reliable  discrimination. 

•  Discrimination  techniques  based  on  L„  group 
velocity  show  separation  between  earthquake 
and  explosion  data  in  some  tectonic  environ¬ 
ments  but  not  others  and  need  further  refine¬ 
ment  and  testing. 


2.2  COMPARISON  OF  REGIONAL  PHASES  FROM  WESTERN  UNITED 
STATES  EXPLOSIONS  AND  EARTHQUAKES 


2.2.1  Introduction 

There  has  been  a  great  deal  of  interest  in  the  potential 
use  of  regional  phase  measurements  as  discriminants  between 
earthquakes  and  underground  nuclear  explosions.  Several  studies 
(e.g.,  Blandford  and  Hartenberger ,  1978;  Blandford,  et  al . , 

1980;  Pomeroy  and  Chen,  1980;  Nuttli,  1981;  Gupta  and  Burnetti, 
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1981)  have  focused  on  the  propagation  characteristics  and  rela¬ 
tive  excitation  of  the  Lg  and  Pg  phases  mainly  because,  in  many 
tectonic  environments,  they  are  the  strongest  regional  phases 
and  therefore  the  most  likely  to  produce  measurable  signals  for 
even  small  events.  The  potential  of  regional  phases,  and  in 
particular  the  Lg  and  Pg  phases,  has  remained  indeterminant 
largely  because  a  definitive  data  base  with  earthquake  and  ex¬ 
plosion  sources  in  close  proximity  has  not  been  assembled  and 
theory  has  not  developed  to  the  point  of  permitting  extrapola¬ 
tion  of  experience  between  different  source  and  propagation  en¬ 
vironments.  In  the  following,  we  focus  on  the  development  of 
a  more  definitive  data  base  to  better  establish  excitation 
differences  for  regional  phases  from  earthquakes  and  under¬ 
ground  nuclear  explosions. 

2.2.2  Analysis  of  NTS  Data 

Blandford  and  Klouda  (1980)  recently  reported  on  the 
Lg,  Pg  and  Pn  phases  from  a  large  group  of  NTS  underground 
nuclear  explosions  recorded  at  the  Tonto  Forest  Observatory 
(TFO)  about  500  km  away.  Using  rather  simple  peak-amplitude, 
time  domain  measurements  for  velocity  windows  corresponding  to 
the  appropriate  regional  phases,  they  developed  relations 
between  Lg,  Pg  and  Pn  amplitudes  (vertical  component)  for  the 
NTS  explosions.  Their  results  for  Pg  and  Lg  are  shown  in 
Figure  14  where  the  slope  of  the  linear  regression  was  con¬ 
strained  to  one  and  the  dashed  lines  represent  2a  error  bounds. 

The  equation  of  the  line  is  log  Ap  =  log  Al_  +  0.1. 

3  y 

In  the  current  program,  Bennett  and  Murphy  (1980)  studied 
a  sample  of  21  earthquakes  covering  a  body-wave  magnitude  range 
from  3.3  to  4.6  which  occurred  within  about  one  degree  of  NTS 
(cf.  Figure  15)  and  were  recorded  at  TFO.  To  avoid  the  ambi¬ 
guity  of  collapse  events  or  unannounced  tests,  we  excluded  all 
events  located  on  the  test  site  except  for  two  events  known  to 
be  earthquakes  near  Massachusetts  Mountain  (m^  ~  4.3)  and  near 

23 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


Figure  15.  Earthquake  sample  from  the  vicinity  of  NTS  with 
regional  phases  recorded  at  TFO.  Circles  indi¬ 
cate  epicenters;  triangle  is  the  station  loca¬ 
tion. 
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Ranger  Mountain  (m^  «  4.2).  We  followed  a  procedure  *  that 
used  by  Blandford  and  Klouda  (1980)  and  measured  peak  amplitudes 
in  time  windows  corresponding  to  the  propagation  velocities  of 
the  expected  regional  phases.  The  necessity  of  such  a  procedure 
is  illustrated  in  Figure  16  which  shows  a  typical  time  history. 
The  maximum  amplitudes  are  somewhat  lacking  in  character  and  in 
many  cases  are  nearly  indistinguishable  from  the  surrounding 
wave  trains.  Figure  17  summarizes  our  results  for  P  and  L 

g  g 

phases  from  the  earthquakes  near  NTS.  The  linear  regression 
was  again  constrained  to  have  a  slope  of  one  resulting  in  the 
relation: 

lo9  APg  -  log  ALg  -  0.1 

for  the  earthquakes.  As  can  be  seen  in  the  figure,  most  of  the 
earthquake  data  points  fall  below  the  mean  explosion  line.  How¬ 
ever,  the  scatter  in  the  earthquake  data  is  so  large  that  the 
2 a  error  bounds  envelop  the  mean  explosion  line  and,  in  fact, 
most  of  the  explosion  data.  For  the  two  earthquakes  at  NTS, 
the  Massachusetts  Mountain  earthquake  data  point  falls  well 
below  the  mean  earthquake  line  and  hence  discriminates  very  well 
from  the  explosion  data.  On  the  other  hand,  the  Ranger  Mountain 
earthquake  data  point  falls  very  close  to  the  mean  explosion 
line  and  would  not  discriminate. 

Similar  measurements  of  the  P  and  L  phases  for  the 

n  g 

earthquakes  and  a  sample  of, the  explosions  in  the  magnitude 
range  from  3.5  to  4.8  is  shown  in  Figure  18.  The  data  appear 
to  be  intermingled  and  offer  little  promise  for  their  use  in 
discrimination. 
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Figure  17.  Comparison  of  Pg  and  Lg  amplitudes  recorded  at  TFO 
from  earthquakes  near  NTS.  Triangles  indicate  data 
points  and  the  solid  line  is  a  regression  line 
through  the  data. 
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2.2.3  Factors  Affecting  the  Discriminants 


Among  the  factors  which  have  been  thought  to  affect  the 
amplitudes  of  regional  phases,  and  hence  to  be  a  source  of  the 
observed  scatter,  are  source  depth,  regional  propagation  path 
and  recording  site  response.  We  will  say  more  about  the  source- 
depth  effect  in  Section  2.3.  In  our  study  of  the  NTS  earthquake 
data  we  attempted  to  minimize  any  differences  in  regional  propa¬ 
gation  path  and  recording  site  by  using  only  a  single  station 
and  limiting  the  earthquakes  to  a  relatively  narrow  band  of 
azimuths  relative  to  the  recording  site  at  TFO .  In  a  second 
part  of  the  study  we  attempted  to  investigate  these  factors 
using  earthquakes  mainly  in  California  (cf.  Figure  19)  recorded 
by  short-period  instruments  in  three  different  areas  of  NTS. 

We  made  measurements  of  the  P  and  L  amplitudes  like 

g  g 

those  which  we  made  above  for  the  TFO  recordings.  The  results 
are  summarized  in  Figure  20.  Different  symbols  correspond  to 
different  recording  sites  on  NTS.  For  each  earthquake  an  ob¬ 
servation  was  generally  available  for  Climax  Stock  and  for 
either  Pahute  Mesa  or  Yucca  Flats.  The  Climax  Stock  observa¬ 
tions  are  confined  mainly  to  lower  values  while  Pahute  Mesa 
and  Yucca  Flats  data  cluster  at  higher  values  suggesting  that 
the  Pahute  Mesa  and  Yucca  Flats  sites  tend  to  amplify  the  P 

g 

and  L  phases  relative  to  Climax  Stock.  A  similar  result  was 

g 

reported  by  Barker,  et  al.  (1980)  using  the  same  data  base. 

It  is  also  interesting  to  note  that  the  trend  in  the  data  re¬ 
mains  close  to  a  slope  of  one.  It  is  unclear  whether  differ¬ 
ences  in  source  excitation  or  attenuation  are  responsible  for 
the  slightly  lower  intercept  of  'hese  data  relative  to  the  TFO 
mean  lines. 


2.2.4  Summary 

Analysis  of  regional  phases  from  earthquakes  and  explo¬ 
sions  in  the  vicinity  of  NTS  and  from  a  supplemental  set  of 
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Figure  20.  Comparison  of  P_  and  Lg  amplitudes  from  California 
earthquakes  recorded  at  different  locations  on  NTS. 
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California  earthquakes  recorded  at  NTS  has  led  us  to  conclude 
that: 

•  The  P  to  L  amplitude  ratios  for  earthquakes 
in  the  vicinity  of  NTS,  based  on  simple  peak 
motion  measurements,  average  less  than  one. 

•  This  observation  contrasts  with  similar  ampli¬ 
tude  ratios  measured  for  NTS  underground 
nuclear  tests  which  average  more  than  one. 

•  The  earthquake  amplitude  measurements  show 
relatively  large  scatter  which  envelopes  the 
explosion  values  and  places  in  question  the 
reliability  of  the  discriminant. 

•  The  relation  between  Pg  and  Lg  amplitudes  from 
earthquakes  does  not  appear  to  be  particularly 
sensitive  to  detailed  differences  in  site  response 
but  may  show  some  dependence  on  transmission  path 
and  focal  depth  (i.e.,  see  following  section). 

•  The  indistinct  character  of  many  of  the  observed 
regional  phases  suggests  that  measurements  other 
than  simple  peak  amplitude  observations  may  better 
represent  signal  strength  in  regional  phase  wave 
trains  and  may  produce  a  better  discriminant. 


2.3  DEPENDENCE  OF  Lg  EXCITATION  ON  SOURCE  DEPTH 
2.3.1  Introduction 

As  was  noted  in  Section  2.2.3  above,  one  of  the  factors 
which  can  be  expected  to  have  an  effect  on  the  amplitudes  of 
regional  phases  is  source  depth.  This  is  especially  true  with 
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regard  to  Lg  which  is  now  generally  recognized  to  correspond  to 
the  superposition  of  higher  mode  Love  and  Rayleigh  surface  waves 
which  are  known  to  have  a  pronounced  dependence  on  focal  depth. 

In  particular,  surface  wave  excitation  is  generally  expected  to 
decrease  with  increasing  source  depth  and,  since  the  amplitudes 
of  regional  P  phases  such  as  Pg  and  Pn  are  known  to  be  relatively 
insensitive  to  focal  depth,  P/Lg  ratios  might  be  expected  to  in¬ 
crease  with  the  depth  of  the  source.  This  suggests  that  deeper 
earthquakes  may  have  explosion-like  P/Lg  ratios,  perhaps  provid¬ 
ing  a  partial  explanation  for  the  overlap  in  the  explosion  and 
earthquake  data  sets  such  as  that  shown  in  Figures  17  and  18. 

We  have  been  investigating  this  issue  using  both  the  results  of 
theoretical  modeling  and  analyses  of  regional  phase  data  re¬ 
corded  at  NTS  from  a  number  of  southern  California  earthquakes 
covering  a  representative  range  in  focal  depths . 

2.3.2  Summary  of  Theoretical  Modeling  Results 

The  Lg  phase  is  a  high  frequency  wavetrain  which  arrives 
with  a  group  velocity  near  3.5  km/sec.  It  is  observed  on  all 
three  components  of  motion  for  both  explosions  and  earthquakes 
and  is  often  the  dominant  arrival  on  short  period  recordings  in 
the  regional  distance  range.  As  a  result  of  its  prominence,  it 
has  been  the  subject  of  numerous  empirical  investigations.  How 
ever,  attempts  to  synthesize  this  phase  theoretically  date  back 
only  a  few  years  and,  consequently,  the  analytical  definition 
of  the  dependence  of  Lg  excitation  on  variables  such  as  source 
depth  and  crustal  structure  is  the  subject  of  ongoing  research. 
Bache,  et  al.  (1980)  recently  completed  a  study  in  which  they 
attempted  to  theoretically  synthesize  vertical  component  Lg 
waveforms  using  superposition  of  Rayleigh  modes  in  a  plane¬ 
layered  crustal  model.  The  applicability  of  this  proposed  syn¬ 
thesis  technique  was  assessed  through  comparisons  with  regional 
data  measured  in  the  eastern  United  States  from  the  SALMON 


34 

SYSTEMS.  SCIENCE  AND 


'JM  IPP|P» jl 


explosion.  Although  these  detailed  comparisons  are  not  entirely 
satisfactory  at  the  present  time,  they  do  suggest  that  the 
theoretical  model  can  explain  the  major  features  of  the  ob¬ 
served  Lg  phase  and  thus  should  be  useful  for  the  parametric 
evaluation  of  the  dependence  of  Lg  excitation  on  the  model 
parameters . 

In  particular,  with  regard  to  the  dependence  of  Lg  ex¬ 
citation  on  focal  depth,  Bache,  et  al.  (1980)  have  synthesized 
Lg  wavetrains  for  point  double  couple  earthquake  sources  at  ten 
different  depths  (i.e.  1,  2,  4,  5,  6,  8,  11,  14,  15  and  25  km) 
in  an  approximate  eastern  United  States  crustal  model.  For 

this  application,  the  seismic  moment  of  the  double  couple  was 
25 

held  fixed  at  10  dyne-cm  and  the  source  time  history  was 
taken  to  be  a  step  function  so  that  the  source  function  has  no 
characteristic  corner  frequency.  Seismograms  were  computed 
for  this  source  oriented  to  represent  strike-slip,  vertical 
dip-slip  and  45°  dip-slip  faulting  (i.e.  the  solutions  for  any 
other  fault  orientation  can  be  obtained  by  linearly  superposing 
solutions  for  these  three  basic  orientations) .  It  was  found 
that  there  is  little  systematic  dependence  of  the  waveforms  on 
depth,  with  the  possible  exception  that  the  duration  of  the 
vertical  dip-slip  synthetics  appears  to  decrease  somewhat  with 
increasing  depth.  However,  the  maximum  amplitude  of  Lg  does 
show  a  pronounced  depth  dependence.  Figure  21  shows  the 
theoretical  Lg  amplitudes  (as  seen  on  the  output  of  a  WWSSN 
short  period  instrument  at  a  range  of  1000  km)  as  a  function 
of  source  depth  for  the  three  source  orientations.  The  source 
locations  are  shown  with  respect  to  the  adopted  crustal  model 
(shear  wave  velocity  S,  attenuation  parameter  Q)  and,  from  this 
display,  it  appears  that  there  is  no  evidence  that  the  Lg 
amplitude  depends  on  the  layer  in  which  the  source  is  located. 
Thus,  to  first  order,  the  depth  dependence  should  be  independent 
of  the  details  of  the  crustal  model.  It  can  be  seen  that  for 
both  the  strike-slip  and  45°  dip-slip  sources  the  amplitude  of 
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Lg  decreases  with  increasing  source  depth,  as  expected.  Over 
the  depth  range  from  1  to  25  km,  the  Lg  amplitude  is  predicted 
to  decrease  by  a  factor  of  four  to  five  for  these  two  source 
orientations.  However,  the  depth  dependence  of  the  Lg  phase 
radiated  from  the  vertical  dip-slip  source  shows  a  more  com¬ 
plicated  pattern  with  the  predicted  amplitude  actually  in¬ 
creasing  with  depth  over  some  ranges. 

Thus,  the  Lg  depth  dependence  predicted  by  available 
theoretical  models  depends  on  source  orientation  and  it  follows 
that  although  a  decrease  in  L amplitude  with  depth  can  be 
expected  in  most  cases,  variations  in  focal  mechanism  may  in¬ 
troduce  considerable  scatter. 

2.3.3  Analysis  of  Southern  California  Earthquake  Data 

In  parallel  with  the  ongoing  theoretical  studies,  we 

have  also  been  analyzing  recorded  earthquake  data  in  an  attempt 

to  empirically  define  the  effects  of  focal  depth  on  the  P/L  type 

y 

regional  discriminants.  For  this  purpose  we  have  selected  a 
sample  of  16  southern  California  earthquakes  with  focal  depths 
(NEIS)  ranging  from  1  to  16  km.  All  of  these  events  were  re¬ 
corded  at  the  same  temporary  station  located  at  Climax  Stock  on 
NTS.  Thus,  the  propagation  path  and  local  receiver  site  effects 
should  be  approximately  constant  and  should  not  mask  any  systema¬ 
tic  effects  of  focal  depth.  The  amplitudes  of  the  various 
regional  phases  recorded  at  the  Climax  Stock  station  from  these 
earthquakes  have  been  read  following  the  procedure  described 
previously  in  Section  2.2.2  and  Figure  22  shows  the  resulting  Pg 
to  Lg  amplitude  ratios  plotted  as  a  function  of  focal  depth.  It 
can  be  seen  that  although  the  data  show  considerable  scatter, 
there  does  appear  to  be  a  tendency  for  the  Pg/Lg  ratio  to  in¬ 
crease  with  increasing  focal  depth,  in  qualitative  agreement 
with  the  theoretical  prediction  for  strike-slip  faulting.  More 
specifically,  for  the  cluster  of  events  near  6  km  focal  depth 
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Figure  22.  Dependence  of  Pg/Lg  amplitude  ratios  on  focal 
depths  for  Southern  California  earthquakes 
recorded  at  NTS. 
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the  amplitude  ratios  average  0.4  to  0.5  while  the  data  at  15 
to  16  km  focal  depths  show  amplitude  ratios  of  near  0.8. 

2.3.4  Summary 

Analyses  of  both  the  results  of  theoretical  modeling  and 
a  limited  sample  of  empirical  data  suggest  that  the  excitation 
of  the  Lg  phase  generally  decreases  with  increasing  focal  depth. 

In  particular,  the  theoretical  simulations  and  observed  data 
both  indicate  about  a  factor  of  two  reduction  in  Lg  peak  amplitude 
level  for  strike-slip  sources  over  the  depth  range  from  about  5 
to  15  km.  This  effect  is  larger  than  the  average  difference 
between  earthquake  and  explosion  Lg  amplitude  levels  and  tends 
to  make  deeper  earthquakes  look  explosion -like  with  respect  to 
the  P/L^  discriminants.  These  preliminary  conclusions  are 
currently  being  tested  against  a  larger  sample  of  regional 
data  recorded  from  Eurasian  events. 


2.4  DEVELOPMENT  AND  TESTING  OF  SPECTRAL  DISCRIMINANTS 

I  2.4.1  Introduction 

In  Sections  2.1  to  2.3,  we  emphasized  the  development 
of  discrimination  techniques  based  on  time  domain  measurements 
(amplitude  and  period)  of  regional  seismic  phases.  In  this 

*  section  we  focus  on  spectral  methods.  In  the  following  dis¬ 
cussion  we  consider  some  preliminary  results  of  the  variable 
frequency  magnitude  (VFM)  approach  (Archambeau,  et  al. ,  1974) 
from  a  study  of  the  SALMON  explosion  and  Alabama  earthquake 

*  pair.  This  work  is  designed  to  test  the  VFM  discrimination 
approach  for  a  closely  located  explosion-earthquake  pair  in 
the  eastern  United  States,  where  relative  attenuation  effects 
are  reduced  and  common  stations  are  available.  We  conclude 

1  this  section  with  a  description  of  a  signal  analysis  approach 


39 


SYSTEMS.  SCIENCE  AN O  SOFTWARE 


which  is  expected  to  provide  more  accurate  spectral  measure¬ 
ments  of  discrimination  parameters  for  isolated  regional  phases. 


2.4.2  Preliminary  Results  for  the  SALMON  Explosion  and 

Alabama  Earthquake  Events 

The  SALMON  and  Alabama  events  were  seismically  recorded 
at  many  common  LRSM  stations.  However,  digital  data  for  the 
Alabama  earthquake  earthquake  are  available  at  S3  for  only  some 
stations  common  to  SALMON.  These  data  lie  along  a  northeasterly 
line  from  the  explosions  and  earthquake  epicenters  (Figure  23) . 
Since  the  objective  of  the  study  is  to  look  at  detailed  source 
spectral  differences  we  need  to  account  for  path  and  receiver 
effects.  This  can  be  accomplished  either  experimentally  by 
selecting  common  paths  and  receiver  combinations  or  theoreti¬ 
cally  treat  the  problem  using  detailed  knowledge  of  crust  and 
upper  mantle  velocity  and  Q  structures.  The  latter  information 
is  known  generally  but  not  in  sufficient  detail  to  adequately 
predict  the  desired  effects.  However,  with  the  experimental 
approach  path  effects  can  be  substantially  negated  by  comparing 
data  along  common  paths  and  equivalent  distances.  That  is,  the 
spectral  difference  observed  should  be  due  to  only  the  source 
if  path  and  receiver  are  the  same.  However,  experimental  real¬ 
ities  prevent  such  ideal  comparisons.  For  example,  digital 
data  are  available  at  EUAL  for  the  Alabama  earthquake  but  not 
for  SALMON .  Digital  data  are  available  at  JELA  for  SALMON  but 
not  for  the  Alabama  earthquake.  None  of  the  paths  are  the  same 
but  both  stations  are  located  in  the  same  general  coastal  crust 
and  upper  mantle  region  (Figure  23).  Thus,  for  short  distances 
the  best  alternative  is  to  directly  compare  JELA  (SALMON)  with 
EUAL  (Alabama  earthquake),  at  distances  of  2.19°  and  2.82°, 
respectively . 

In  Figures  24  and  25  we  show  arrival  times  of  various 
phases  generated  from  dispersive  filtering  of  signals  at  JELA 
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Figure  23.  Station  location,  signal  reception,  structural  regions,  and 
station-source  profiles. 


i  • 


Figure  24.  JELA,  SALMON  seismogram  showing  various  phase 

arrival  times  generated  by  dispersion  filtering 
using  only  "one  pass"  detection  (A  *  243  km, 
STA-EPI  azimuth  »  106°). 
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Figure  25.  EUAL,  Alabama  earthquake  seismogram  showing 

various  phase  arrival  times  generated  by  dis¬ 
persion  filtering  using  "one  pass"  detection 
only  (A  =  314  km,  STA-EPI  azimuth  =  44°). 
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and  EUAL.  The  signals  appear  quite  different  in  shape  and  size 
but  not  in  frequency  content.  This  is  evident  from  narrow- 
band  filter  amplitude  estimates  shown  in  Figures  26  and  27. 
Figure  26  compares  the  maximum  ?  pulse  amplitude  spectra  for 
the  two  events  and  Figure  27  shows  the  maximum  P^  pulse  spectra. 

To  scale  the  spectra  at  a  common  distance,  the  JELA/ 

SALMON  P  spectrum  is  divided  by  a  factor  of  two  based  on  an 
n  -3 

attenuation  rate  of  A  ^  r  for  the  western  United  States  (WUS) 

(Blandford,  et  al.,  1980).  The  attenuation  of  P  in  the  eastern 

—  —  n  -2  5 

United  States  is  probably  somewhat  less,  perhaps  A  'v-  r 

(Blandford,  et  al . ,  1980)  but  a  factor  of  two  is  adequate  for 

the  spectral  comparison  particularly  since  small  differences 

in  attenuation  (t*)  which  may  occur  throughout  the  coastal 

region  structure  will  not  significantly  change  the  character 

of  the  spectral  differences.  The  Pg  spectra  for  SALMON  is 

divided  by  a  factor  of  1.5.  This  factor  is  based  on  the 

attenuation  rate  for  Pg  derived  by  Baker,  1970;  also  for  the 

WES. 

We  observe  that  the  SALMON  Pn  spectrum  is  greater  in 

value  and  peaks  at  a  higher  frequency  than  the  spectrum  for 

the  earthquake  and  this  is  true  even  if  the  spectra  is  averaged 

or  smoothed.  The  P  spectral  estimates  (for  the  P„  maximum 

g  g 

pulse)  from  the  two  events  are  quite  similar.  This  similarity 

in  the  maximum  Pg  pulse  may  be  due  more  to  the  bandpass  char¬ 
acter  of  the  structure  transfer  function  than  to  similarity  in 
the  two  source  spectra.  The  former  observation,  however,  sug¬ 
gests  that  the  variable  frequency  magnitude  discriminant 
applied  to  PR  may  be  effective  for  regional  events.  On  the 
other  hand,  a  spectral  discriminant  based  on  the  maximum  Pg 
pulse  may  not  be  useful.  Many  more  station-event  pairs  will 
be  required  to  verify  these  preliminary  observations  however. 

Figures  28  and  29  show  the  seismograms  and  various 
arrival  times  obtained  irom  dispersion  filtering  applied  to 
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Figure  28.  BLWV,  SALMON  seismogram  showing  various  phase 
arrival  times  derived  by  first  pass  dispersion 
filtering  (A  =  1058  km,  STA-EPI  azimuth  =*  228°)  . 
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Figure  29.  BR-PA,  Alabama  earthquake  seismogram  showing 
various  phase  arrival  times  derived  by  first 
pass  dispersion  filtering  (1  =  818  km,  STA-EPI 
azimuth  =*  228°)  . 
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the  two  stations  BLWV  and  BRPA,  located  at  distances  of  1058 
and  818  kilometers  from  SALMON  and  the  earthquake,  respectively. 
The  signals  arrive  at  both  stations  from  the  same  direction 
(i.e.,  STA-EPI  azimuth  =  228°).  This  event-station  pair  was 
selected  for  nearly  equal  distances  and  azimuths,  again  to 
minimize  path  effects.  Granted  the  receiver  and  path  transfer 
functions  may  differ  (Figure  23)  and  a  better  selection  would 
be  SALMON/BLWV  and  earthquake/DHNY  as  far  as  distance  and  path 
are  concerned,  but  again  experimental  realities  interfer. 

That  is  the  signal-to-noise  (S/N)  ratio  is  too  low  for  the 
earthquake/DHNY  record  and  as  will  be  shown  the  S/N  ratio  is 
also  low  for  BRPA. 

The  Pn  phase  at  BRPA  for  the  earthquake  (Figure  29)  is 
obscured  by  noise;  however  a  larger  signal  is  present  for  the 
P  phase.  In  Figure  30,  we  compare  the  P  peak  amplitude 
spectra  for  the  two  events.  The  spectrum  for  BLWV/ SALMON  is 
divided  by  1.5  to  adjust  for  the  existing  distance  difference 
(Baker,  1970) .  The  BRPA  Pg  earthquake  spectrum  peaks  at  a 
slightly  lower  frequency  than  that  at  BLWV  for  the  explosion, 
and  generally  shows  larger  amplitudes  over  the  band  of  fre¬ 
quencies  analyzed.  However,  in  view  of  the  rather  low  signal- 
to-noise  ratio,  it  is  probably  not  realistic  to  place  much 
weight  on  these  rather  small  differences.  Rather,  it  seems 
most  appropriate  to  investigate  the  possible  Pg  phase  dif¬ 
ferences  between  the  source  types  over  more  of  the  wavetrain 
(that  is  using  more  than  just  the  maximum  pulse  within  the 
wavetrain)  and  over  a  wider  bandwidth,  in  particular  at  lower 
frequencies.  These  investigations  are  in  progress,  using  more 
sophisticated  quasi-harmonic  decomposition  methods. 

In  addition  to  the  P„  and  P„  phases,  the  L„  "phase"  is 

n  g  g 

a  possible  discriminatory  arrival.  Both  P_  and  L  can  be 

9  9 

described  as  the  superposition  of  a  number  of  higher  order  P- 
SV  or  Rayleigh  type  modes  (Press  and  Ewing,  1952  and  Knopoff, 
et  al . ,  1973).  Alternatively,  both  can  be  described  in  terms 
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Figure  30.  NBF  spectra  for  the  maximum  Pg  phase  (U  <  6.0  km/sec)  at  BRPA  and 
BLWV  for  the  Alabama  earthquake  and  the  SALMON  explosion. 


of  a  large  number  of  multiple  reflected  body  phases  controlled 
primarily  by  P  or  S  velocities,  respectively  (Archambeau,  1981; 
Helmberger  and  Enger,  1980) . 

Figure  31  shows  the  arrival  times  of  pulses  within  an 
Lg  group  velocity  window  of  3.6  to  2.9  for  the  Alabama  earth¬ 
quake,  recorded  at  BRPA.  The  spectrum  of  the  maximum  pulse 
in  the  Lg  wavetrain  is  plotted  in  Figure  32 ,  along  with  the 
maximum  Pg  pulse  spectrum  and  the  noise  spectrum  estimate  for 
this  event.  The  Lg  spectrum  peaks  at  a  slightly  lower  fre¬ 
quency  than  Pg  and  is  richer  in  low  frequency  energy.  On  the 
other  hand,  from  preliminary  time  domain  observations,  the 
explosion  event  data  appears  to  show  Lg  less  than  or  comparable 
to  Pg  at  most  stations.  This  suggests  that  the  relative 

excitation  of  the  maximum  P  to  L  may  serve  as  a  useful  dis- 

g  g 

criminant.  We  currently  are  applying  advanced  signal  isola¬ 
tion  and  spectral  estimation  methods  (see  Section  2.4.3  for 
example)  to  systematically  verify  these  very  preliminary 
observations. 

The  effects  of  receiver  and  source  structure  as  well  as 
crust-mantle  characteristics  on  the  ground  motion  at  regional 
distances  have  been  discussed  by  many  authors.  Additional, 
systematic  characterization  of  these  effects  is  indeed  an 
important  part  of  this  study.  However,  because  of  the 
complex  variations  and  indistinct  character  of  many  of  the 
observed  regional  phases,  measurements  of  other  than  simple 
peak  amplitude  spectral  estimates  will  almost  certainly  be 
required  for  the  development  of  robust  regional  discriminants. 
In  this  regard,  energy  estimates  derived  for  isolated  pulse 
arrivals  and  for  the  entire  wavetrain  are,  and  will  be,  in¬ 
vestigated. 

By  way  of  example,  simple  integration  of  the  wavetrain 
spectra  over  a  frequency  band  for  a  specified  time  interval 
would  tend  to  average  out  the  effects  of  variations  in  propa¬ 
gation  over  different  paths.  For  example,  Figures  33  and  34 
show  Lg  signals  from  two  explosions,  one  in  the  EUS  (SALMON) 
and  the  other  in  WUS  (DIABLO  HAWK)  recorded  at  EUAL  and 
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31.  BRPA,  Alabama  earthquake  seismogram  showing 

various  phase  arrival  times  generated  by  dis¬ 
persion  filtering  (A  *  818  km,  STA-EPI  azimuth  * 
228°) . 
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Minersville,  DIABLO  HAWK  seismogram  for  a  group  velocity  window  of 
3.5  to  2.4  km/sec  (A  =  303  km,  STA-EPI  azimuth  =  250°). 


Minersville,  Utah,  respectively.  For  the  EUS  (Figure  33) ,  the 
1.^  energy  is  relatively  greater  in  the  later  portion  of  the 
wavetrain  than  for  the  early  portion.  The  opposite  is  true 
for  the  WUS  recording  (Figure  34) .  The  maximum-entropy  moving- 
window  spectra  confirms  these  observations  (Figures  35  and  36) . 
The  maximum-entropy  moving  window  spectra  is  similar  to  a  con¬ 
ventional  sonogram  analysis  described  by  Cohen  (1969) .  We  ex¬ 
pect,  however,  that  proper  summations  over  time  would  be 
independent  of  time  shifts  in  the  energy  content  due  to  propa¬ 
gation  effects,  and  provide  a  stable  estimate  of  the  total  re¬ 
corded  energy  in  the  P  and  L  wavetrains.  In  this  case,  we 

g  g 

expect  that  the  relative  energy  content  of  P_  to  L  will  pro- 

9  9 

vide  a  good  event  discrimination  variable.  Our  current  and 
proposed  work  is  directed  to  a  systematic  investigation  of 
this  and  similar  possibilities,  as  suggested  by  this  preliminary 
work. 


A  more  sophisticated  form  of  the  signal  detection/ 
analysis  procedure  has  been  developed  (Archambeau,  et  al. , 
1981a)  and  is  now  incorporated  in  the  MARS  (Multiple  Arrival 
Recognition  System)  program.  We  generally  refer  to  the  pro¬ 
cedures  involved  in  this  time  series  analysis  approach  as 
multiple-pass  dispersion  filtering.  The  entire  approach  is 
based  on  quasi-harmonic  decomposition  of  the  time  series, 
using  narrow  band  filtering  techniques,  and  is  described 
quantitatively  by  Archambeau,  et  al . ,  1981b.  The  approach 
has  great  flexibility  in  that,  when  three  component  seismic 
data  is  available,  both  polarization  and  wave  number  filtering 
can  be  carried  out  jointly  with  the  dispersion  filtering  to 
achieve  a  high  degree  of  selectivity  for  signal  detection 
and  a  large  amount  of  signal  information  for  all  the  signals 
detected.  For  regional  data,  where  a  large  number  of  seismic 
phases  arrive  close  to  each  other  in  time,  it  is  not  only 
desirable,  but  quite  necessary,  to  isolate  individual  phases 
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recorded  at  EUAL  from  the  SALMON  explosion  (contour  interval  =  5  dB 
242  km,  STA-EPI  azimuth  =  222°). 


Figure  36.  Contoured  maximum-entropy  moving  window  spectral  estimates  of  the  Lq  phase 
recorded  at  Minersville,  Utah  from  the  DIABLO  HAWK,  NTS  explosion  (contour 
interval  =  5  dB,  A  =  303  km,  STA-EPI  azimuth  =  250°). 


and  obtain  good  (uncontaminated)  estimates  of  the  spectral 
content  of  each,  their  individual  arrival  times  and  sense 
of  first  motion,  and,  ideally,  also  their  polarization,  dis¬ 
persion  and  wave  number  vector  orientations  as  functions  of 
frequency. 

The  new  computer  program  developed  has  the  capability 
of  performing  the  entire  complement  of  signal  detection/ 
analysis  procedures  mentioned.  In  the  present  framework,  we 
are  employing  it  to:  (1)  generate  the  signal  data  needed  to 
define  discrimination  variables  and  location  parameters  for 
regional  events,  and  (2)  to  investigate  the  robust  character 
of  candidate  discrimination  variables,  using  both  real  and 
synthetic  seismic  data,  as  functions  of  noise  levels,  event 
type  and  magnitude  and  structural  variations.  Ultimately, 
we  also  expect  to  use  this  program  to  automatically  generate 
a  subset  of  required  discrimination  variables  to  be  used  in  a 
multivariate  discrimination  procedure  for  both  regional  and 
teleseismic  events  (Farr’li,  et  al. ,  1981). 

Figures  37  through  39  provide  an  example  of  the  testing 
of  this  program  for  signal  isolation  and  spectral  analysis 
using  regional  synthetic  seismograms.  The  objective  here  is 
to  verify  that  the  filtering  correctly  isolates  regional  type 
signal  pulses  in  the  presence  of  noise  and  obtains  proper 
spectral  estimates  for  each  known  signal  pulse  in  the  wave- 
train.  Figure  37  shows  predicted  ground  motion  as  a  function 
of  distance  from  the  source,  in  the  near-regional  distance 
range.  The  source  was  chosen,  for  this  test,  to  be  an  explo¬ 
sion  at  great  depth  (10  km) ,  so  that  the  synthetic  seismograms 
would  be  relatively  simple.  The  method  used  to  generate  the 
synthetics  is  the  locked  mode  approximation  devised  by  Harvey, 
1981.  The  phases  arriving  coherently  across  this  distance 
range  can  quite  easily  be  identified,  so  that  at  a  given  dis¬ 
tance,  it  is  possible  to  identify  particular  reflected,  re¬ 
fracted  and  direct  phases.  (Phase  identifications  can  also 
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receiver  range 


Time  (sec) 


Figure  37.  Synthetic  near  regional  distance  seismograms 
generated  by  mode  superposition,  for  testing 
of  signal  detection  and  analysis  methods. 
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Figure  39.  Signal  pulse  spectra  as  estimated  by  QHD  filtering 
methods  from  the  synthetic  seismogram  of  Figure  37. 
The  points  are  directly  estimated  from  multiple 
narrow  band  filter  output  while  the  lines  denote 
the  interpolation-extrapolation  of  this  data  over 
the  band  0  to  12  Hz. 
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be  checked  by  theoretical  ray  tracing  through  the  model  veloc¬ 
ity  structure  used.) 

Figure  38  shows  the  results  of  adding  seismic  noise  to 
the  synthetic  seismogram  at  50  km  from  the  source  (top  trace) 
and  then  processing  the  resulting  time  series  through  the 
dispersion  filter  procedure.  In  this  operation,  only  pulse¬ 
like  signals  with  little  or  no  dispersion  (i.e.,  body  waves) 
are  selected  as  signals  of  interest.  The  procedure  involves 
ultra-narrow  band  filtering  and  generation  of  envelope  and 
instantaneous  phase  functions  from  which  the  energy  arrival 
times  (group  arrival  times)  and  the  spectral  amplitude  and  phase 
at  these  times  are  measured  and  associated  with  the  center 
frequencies  of  the  narrow  band  filters.  The  desired  undis¬ 
persed  pulses  are  then  selected  by  a  pattern  search  in  the 
frequency- time-amplitude  parameter  space,  with  those  pulses 
having  the  proper  dispersion  and  preselected  frequency  depen¬ 
dent  signal-to-noise  threshold  level  being  selected.  Pulses 
are  selected,  their  spectra  subtracted  from  the  origin  data 
and  the  resulting  "reduced"  origin  spectrum  is  inverted  to  the 
time  domain.  The  entire  filtering  and  search  process  is  then 
repeated.  This  procedure  is  continued  until  no  new  signal 
pulses  are  found.  (This  iteration  procedure  is  the  origin 
of  the  term  "multiple-pass  dispersion  filtering.")  Pulses 
obtained  in  this  way  may  overlap  in  time,  and  it  is  possible 
to  separate  phases  that  are  time  shifted  by  only  a  fraction 
of  their  apparent  period.  On  the  other  hand,  if  the  pulses 
are  significantly  smaller  than  the  uncertainty  or  resolution 
time  of  the  narrow  band  filter,  then  they  are  automatically 
combined  together  (spectrally  added)  at  the  end  of  the  itera¬ 
tion  process  to  form  a  single  pulse.  The  individual  pulse 
spectra  so  isolated  are  then  inverted  to  the  time  domain,  as 
shown  in  the  lower  six  traces  in  Figure  38.  In  the  figure 
these  may  then  be  summed  to  give  a  single  time  series,  as  in 
the  second  from  the  top  trace  in  the  figure,  which  can  be 
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compared  with  the  original.  As  can  be  seen  from  this  compari¬ 
son,  the  detection  and  isolation  of  +-he  pulses  by  the  filter¬ 
ing  operations  appears  to  be  remarkedly  accurate.  This  can 
also  be  verified  by  comparing  the  filtered  output  results 
with  the  original  synthetic  in  Figure  37.  Close  examination 
shows  that  every  significant  pulse  was  detected  and  properly 
isolated  with  no  false  alarms.  (One  very  weak  arrival  was, 
however,  missed  in  the  detection  procedure,  presumably  due 
to  its  very  low  signal-to-noise  ratio.) 

Figure  39  shows  the  spectral  estimates  obtained  for 
each  of  the  pulses  detected.  That  these  spectral  estimates 
are  very  reasonable  estimates  of  the  true  spectra  is  apparent 
from  the  fact  that  when  inverted  to  the  time  domain,  they  re¬ 
sult  in  pulses  that  very  closely  match  those  of  the  pure  syn¬ 
thetic,  without  noise. 

In  addition  to  the  dispersion  and  spectral  estimates  of 
the  isolated  pulses,  it  is  also  possible  to  obtain  the  "times 
of  arrival"  of  the  phases.  These  are  indicated  by  the  dashed 
lines  in  Figure  38. 

Figures  40  and  41  illustrate  the  application  of  the 
procedure  to  a  small  regionally  recorded  earthquake  from 
Southern  California.  Only  the  first  three  seconds  of  the 
time  series  is  shown  in  Figure  40  and  within  approximately 
a  one  second  time  interval  five  separate  overlapping  signal 
pulses  are  detected.  The  separate  pulse  waveforms  are  shown 
in  the  lower  five  traces,  which  have  been  inverted  to  the 
time  domain  from  the  spectral  estimates  shown  in  Figure  41. 

The  sum  of  these  pulses  (second  trace  in  Figure  40)  is  re- 
raarkedly  similar  to  the  first  second  of  the  event  recording. 

The  timing  lines  shown  are  computed  in  a  manner  designed  to 
indicate  where  the  maximum  of  the  early  arriving  energy  occurs, 
rather  than  at  the  signal  pulse  onset  time. 


Time  (sec) 


Figure  40.  Regionally  recorded  seismogram  (A  ^  100  km)  and 
signal  pulses  determined  by  multiple-pass  QHD 
methods  employing  dispersion  filtering.  The  top 
trace  is  the  original  earthquake  data  while  the 
lower  traces  are  time  series  constructed  from 
the  isolated  pulse  spectra.  The  second  trace 
is  the  superposition  of  the  single  pulse  results. 
The  timing  lines  shown  are  determined  automati¬ 
cally  from  (weighted)  mean  group  arrival  times 
for  each  pulse. 
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Figure  40  lists  the  arrival  tines  of  the  isolated 
pulses  and  it  is  evident  that  they  arrive  within  a  few  tenths 
of  a  second  from  each  other. 

It  is  evident  from  these  examples  that  we  are  in  a 
position  to  rapidly  and  accurately  generate  signal  data, from 
which  event  discrimination  and  yield  or  magnitude  results  may 
be  obtained.  In  particular,  we  clearly  can  measure  frequency 
dependent  magnitudes  over  a  wide  band  of  frequencies  for  any 
and  all  isolated  phases.  We  can  also  automatically  determine 
arrival  times  and  first  motion  data  for  all  of  the  phases  ; 

quite  reliably.  We  can,  therefore,  investigate  spectral  dis¬ 
criminants  for  p  ,  P  .  S  ,  L  ,  R  ,  etc.  and,  as  well,  obtain  ! 

n  g  n  g  g 

complexity  measures,  fault  plane  solutions  and  location  and  j 

I 

source  depth  estimates  using  all  the  phase  information  in-  i 

eluding  S-P,  pP,  PP,  etc.  In  this  regard,  we  plan  to  form 
candidate  discriminants  and  systematically  test  them  against  ' 

both  real  and  synthetic  regional  seismic  data.  j 
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III.  SUMMARY  OF  CURRENT  STATUS  AND  FUTURE  STUDIES 

3.1  REGIONAL  PHASE  CHARACTERIZATION 

3.1.1  Introduction 

As  we  described  in  previous  sections  of  this  report, 
development  of  a  reliable  discriminant  using  regional  phase 
amplitude  measurements  has  been  hindered  by  the  substantial 
scatter  in  the  measurements  and  the  lack  of  procedures  for 
extrapolating  experience  from  a  tested  area  to  a  tectonically 
different  untested  area.  In  the  future,  we  will  investigate 
factors  which  might  be  responsible  for  the  observed  scatter, 
in  the  hope  of  developing  corrections  to  the  data  which  will 
permit  a  reduction  in  the  scatter  and  hence  result  in  a  more 
reliable  discriminant.  Those  factors  which  appear  to  contri¬ 
bute  to  the  scatter  include  path  attenuation,  earthquake  focal 
depth  differences,  earthquake  focal  mechanisms,  source  geo¬ 
logic  environment  and  possible  recording  site  conditions.  It 
is  our  intention  to  analyze  these  factors  using  a  well-con¬ 
trolled  set  of  observational  data  and  theoretical  studies  of 
the  dependence  of  regional  phase  excitation  and  propagation 
on  these  factors. 

3.1.2  Attenuation  of  Regional  Phases 

Our  plan  is  to  continue  to  develop  a  data  base  for  de¬ 
scribing  the  behavior  of  regional  phase  amplitude  with  distance 
for  transmission  paths  crossing  the  U.S.S.R.  We  intend  to 
categorize  the  available  data  according  to  tectonic  regions. 

For  this  purpose,  distinction  will  be  made  between  paths  cross¬ 
ing  the  platform  and  older  fold  areas  and  those  crossing  younger 
folded ^areas  which,  on  the  basis  of  our  initial  results,  re¬ 
present  a  clear  difference  in  regional  phase  transmission  prop¬ 
erties.  In  particular,  we  will  investigate  in  greater  detail 
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the  attenuation  of  regional  phases  leaving  the  Soviet  Union 
I  to  the  south  were  stations  important  to  the  identification  of 

small  events  in  known  Soviet  nuclear  testing  areas  may  be 
located . 

Using  the  available  data  we  intend  to  develop  a  more 
1  quantitative  description  of  the  effects  of  attenuation  on  re¬ 

gional  phase  amplitudes.  Such  a  description  would  ideally 
include  identification  of  an  effective  Q  such  as  described  by 
Nuttli  (1980)  for  various  regional  phases  along  paths  from 
specific  areas  of  particular  interest  (e.g.,  known  Soviet 
test  areas)  to  regional  stations  important  to  event  discrimina¬ 
tion  in  those  areas .  Where  it  is  not  possible  to  do  this 
directly  on  the  basis  of  published  data,  it  may  be  possible  s 
to  infer  more  quantitative  Q  estimates  from  regions  of  similar 
tectonics  where  the  attenuation  of  regional  phases  is  known. 

A  firm  basis  for  inferring  attenuation  properties  could  be 
established  from  additional  theoretical  analysis  of  regional 
phase  transmission. 

3.1.3  Variations  in  the  Source 

The  principal  factors  to  be  considered  here  appear  to 
be  source  depth  differences  and  effects  of  focal  mechanism 
on  radiation  pattern.  As  described  in  Section  2.3,  the  focal 
depth  is  expected  to  have  significant  effects  on  the  excita¬ 
tion  of  certain  regional  phases.  This  is  demonstrated  for 
the  Lg  phase  theoretically  by  the  level  of  higher  mode  sur¬ 
face  wave  excitation  from  simple  double  couple  sources  at 
various  depths  of  burial.  We  also  noted  a  possible  depth 

dependence  in  observations  of  the  amplitude  ratios  of  the  L 
i  5 

and  Pg  phases  for  events  in  Southern  Calif orna.  We  are  cur¬ 
rently  in  the  process  of  developing  a  data  set  covering  a  much 
larger  range  in  focal  depths  using  earthquakes  in  the  Pamir- 
Hindu  Kush  area  recorded  at  regional  distances.  It  is  our 


69 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


» 


» 


t 


» 


» 


intention  to  use  measurements  of  regional  phase  amplitudes 
from  these  data  to  better  characterize  the  dependence  of  ex¬ 
citation  of  various  regional  phases  on  source  depth.  For  the 
Lg  phase,  we  will  compare  the  observations  with  theoretical 
computations  of  regional  phase  excitation  based  on  the  higher 
mode  surface  wave  representation. 

Finally,  we  will  investigate  the  possible  effect  of 
variations  in  source  radiation  pattern  in  producing  the  ob¬ 
served  scatter.  Clearly,  for  earthquakes  from  a  restricted 
source  region  we  would  expect  this  effect  to  be  minor.  However, 
in  many  of  the  studies  using  Soviet  data,  there  is  a  large 
separation  between  earthquake  sources  which  make  up  a  common 
data  set.  Wherever  data  permit,  we  intend  to  take  a  closer 
look  at  possible  mechanisms  of  individual  events  which  con¬ 
tributed  to  these  data  sets.  Using  the  higher  mode  excita¬ 
tion  model,  we  will  then  assess  the  significance  of  mechanism 
differences  in  producing  the  observed  scatter. 


3.1.4 


jlication  to  Discrimination 


Using  the  results  of  the  above  studies,  we  plan  to 
analyze  observations  for  a  data  base  including  earthquakes 
and  explosions  in  the  U.S.S.R.  We  will  measure  amplitudes 
of  regional  phases  for  the  events  in  this  data  base.  After 
making  appropriate  corrections  for  factors  believed  to  con¬ 
tribute  to  the  scatter,  we  will  attempt  to  identify  excitation 
differences  between  various  regional  phases  for  earthquakes 
versus  the  explosion  sources.  Finally,  the  reliability  of 
any  discriminant  which  is  proposed  from  the  data  will  be 
assessed  in  terms  of  the  variability  of  the  measurements  and 
the  significance  of  the  correction  factors  which  have  been 
applied . 
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3.2 


SPECTRAL  DISCRIMINANTS  BASED  ON  REGIONAL  PHASES 


One  of  our  future  objectives  is  to  investigate  spectral 
discriminants  using  more  comprehensive  regional  data  sets.  In 
this  regard,  we  will  more  thoroughly  investigate  VFM  type  dis¬ 
criminants,  as  well  as  other  discriminants,  for  regional  data, 
using  large  data  sets  (Caltech  network  data,  LLL  regional  data 
and  AEDS  data  as  described  in  a  separate  section  appended  to 
this  report)  and  synthetically  produced  "data,"  which  can  be 
systematically  analyzed  for  definition  of  discriminants  and 
their  properties  (e.g.,  "robustness")  as  functions  of  super¬ 
posed  noise  levels,  event  type  and  magnitude,  and  source- 
receiver  transmission  characteristics.  Among  the  "VFM"  type 
discriminants  to  be  used,  which  are  essentially  spectral  dis¬ 
criminants  involving  differences  in  spectral  excitation  in  two 
different  narrow  bands  for  the  same  phase,  we  are  considering 
body  wave  magnitude  ("m.  (f)")  discrimination  variables  of  the 
form: 

xk  =  mb(fk)  *  “b^k*  +  bk?  k  *  l’  2'  11  *  k  • 


H  T 

Here  m^Cf^)  and  m^tf^)  denote  magnitudes,  for  a  particular 
body  wave  phase,  at  sets  of  high  (f^)  and  low  (f^)  frequencies, 


and  b^  are  fixed  constants. 


In  general,  a  number  of  sets  of 
high  and  low  frequencies  are  used,  so  that  the  index  k  is  used 
to  label  the  various  combinations  of  high  and  low  frequency 
pairs.  Definition  of  "discrimination  variables"  in  this  way 
is  coupled  to  the  use  of  multivariant  discrimination  methods, 
which  are  described  by  Farrell,  et  al.  (1981). 


Defined  in  this  manner,  we  expect  the  mean  y^tQ)  (or 
median)  of  an  earthquake  "x^  population"  to  be  less  than  the 
mean  y^E)  of  a  similarly  located  explosion  population.  The 
constants  b^  are  used  to  adjust  the  variable  for  regional  bias 
effects. 


71 


SYSTEMS.  SCIENCE  ANO  SOFTWARE 


For  regional  events,  we  are  considering  this  "m^f) 

discriminant  variable"  for  P,  P  ,  P  ,  S  and  S  type  "body 

n  g  n 

wave"  phases. 

In  addition  to  the  body  wave  VFM  discrimination 
variables,  we  also  will  study  similar  surface  wave  type  "VFM 
discriminants."  In  particular: 

=  Ms(f?)  "  Ms(f2}  +  B£;  a  =  1'  2"*;  L  • 

This  very  similar  variable  is  defined  in  a  completely  analogous 
fashion  as  the  "m^ff)  discriminant."  It  is  applied,  for 
regional  events,  to  both  Rayleigh  (R)  and  Love  (L)  type  funda¬ 
mental  mode  surface  waves,  including  Lg  (which  is  composed  of 
higher  order  Rayleigh  and  Love  modes) . 

In  addition  to  the  "VFM  type"  discrimination  variables, 
we  also  consider  the  class  of  what  can  be  termed  the  "mb-Ms 
discriminants."  These  require  mixed  seismic  phase  measure¬ 
ments,  involving  magnitudes  of  two  different  seismic  phases 

such  as  P  and  L  or  P  and  R  at  different  frequencies.  Thus, 
n  g  n 

we  consider  discrimination  variable  sets  of  the  form: 

Z  =  m.  (fj)  -  M  ( f^)  +  C  ;  n  »  1,  2,..,  N 
n  on  s  n  n 

where  the  magnitudes  have  essentially  the  same  definitions  as 
for  the  previously  defined  discriminants.  In  this  case,  how¬ 
ever,  the  magnitudes  maybe  network  averages. 

For  definition  of  Z  for  regional  seismic  data,  we  con- 

**  __  U 

sider  P  .  P  P„  and  S,  S„  for  definition  of  an  m.  (f  ) ,  while 
n  ,  g  _n  -  on 

R,  L  or  L  is  used  for  M.  ( f") .  As  with  the  other  discrimina- 
g  on 

tion  variables,  the  expected  population  means  or  medians  for 
explosions  and  earthquakes  from  similar  regions  is  such  that 

U_  (E)  >  u  (Q)  . 
n  n 
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Finally,  a  number  of  other  discriminants  are  being 

investigated,  in  particular,  complexity  type  discriminants, 

depth  discriminants  involving  time  differences  between  early 

and  later  arriving  phases  in  the  seismogram,  and  time  averaged 

energy  differences  between  the  more  complex  phases,  such  as 

P„  and  L  . 

9  9 
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